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Abstract 



We consider a metric graph Q made of two graphs Qi and Q2 attached at one point. We 
^ ^ derive a formula relating the spectral determinant of the Laplace operator Sg{'y) = det(7— A) 

in terms of the spectral determinants of the two subgraphs. The result is generalized to 
describe the attachment of n graphs. The formulae are also valid for the spectral determinant 



i of the Schrodinger operator det(7 — A + V{x)). 

PACS numbers : 02.70.Hm ; 02. 10. Ox 

J> ' Introduction.— Let us consider a bounded compact domain Vi, part of a manifold. We denote 

. by Spec(— A; Pi) the set of solutions E of —Aip^r) = Eip^r) with ^(r) satisfying given boundary 

conditions at the boundary dVi (Sturm-Liouville problem). Similarly we consider a second 
bounded compact domain 7)2, distinct from Vi and denote Spec(— A;P2) the spectrum of the 

■ Laplace operator in 1)2- Now, if we can glue Pi and D2 by identification of parts of dVi and 
i d'D2 in order to form a unique compact domain D, the question is : can we relate the spectrum 

^ I Spec(— A;2?) to Spec(— A;2?i) and Spec(— A;P2) ? The aim of this article is to discuss this 

question in the particular case of metric graphs when two graphs are attached at one point. 
For that purpose the spectral information is encoded in the spectral determinant of the graph 
^ I Q, formally defined as 5*5(7) = det(7 — A). We first define basic notations and briefly recall 

■ some results on the spectral determinant of metric graphs. We derive the relation between the 
spectral determinant of a graph in terms of the two subgraph determinants, when subgraphs 
are attached by one point, as represented on figure [2jc. The relation is generalized to describe 
attachment of n > 2 graphs (figure [2jd) and to deal with Schrodinger operator. It is interesting 
to point out that our result is reminiscent of the gluing formula for elliptic operators acting on 
a manifold obtained in Ref. |!lj. 

Metric graphs.— Let us consider a collection of V vertices, denoted here by greek letters a, /3 
connected between each others by B bonds, denoted (a/3), {^v)... Each bond is associated with 
two oriented bonds, that we call arcs and denote as a/3, /3a, /xi^, u^.... The topology of the graph 
is characterized by its adjacency (or connectivity) matrix a^^p : a^^p = 1 if (a/3) is a bond and 
dap = otherwise. The coordination number of the vertex a is denoted rua = Ylp^aP- Up 
to now we have built a "combinatorial graph" . If now each bond is identified with an interval 
[0,^0/3] G M, where lap is the length of the bond (a/3), the set of all connected bonds forms a 
"metric graph" (also called a "quantum graph"). 
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A scalar function (p{x) living on a graph Q is defined by B components '■Papixap) where 
Xap £ [0, lap] is the coordinate along the bond {x^^p = corresponds to vertex a and = lap to 
vertex By construction Xap + x^a = lap- Note that components are labelled by arc variables, 
since we must specify the orientation of the axis along which the coordinate is given ; Obviously 
^ap{xap) = ^pa{xf3a) for a Scalar function. The action of the Laplace operator on the scalar 
function along a bond coincides with the one-dimensional Laplace operator {Aip)ap{x) = '^'^pi^)- 
In order to define a self-adjoint operator, one must specify boundary conditions at the vertices. 
The most general conditions have been discussed in Ref. [2] (in general the question of boundary 
conditions is related to the precise nature of the scattering at the vertex O IU O [6]). In the 
present article we consider the simple case of the Laplace operator acting on scalar functions 
that are continuous at the vertices : ^api^ap = 0) = ipa^ P neighbour of a (that gives rua — 1 
equations at the vertex a of coordination number rUa)- Then one must impose another condition 
on derivatives of the function. For continuous boundary condition, the most general condition 
that ensures self-adjointness of Laplace operator is cLap^'api^ap = 0) = ipa with S M. 
The presence of the adjacency matrix in the sum, contraints this latter to run over vertices 
neighbour of a only. The nia equations ensure self adjointness of Laplace operator. Aq, = oo 
corresponds to Dirichlet boundary condition (ipa = 0). The study of the Laplace operator on 
a metric graph appears in several contexts, reviewed in Refs. [Tj [8], like quantum mechanical 
problems : —A(p{x) = E(p{x) could be the Schrodinger equation. In such a case it can be more 
interesting to consider the situation of a graph submitted to a magnetic field, what is achieved 
by replacing the derivative by the covariant derivative, ^ ^ = — ^A{x), where A{x) is 
the vector potential. The boundary condition then reads J2p'^ap(^xV')ap{0) = Xa^a- 




Figure 1: Examples. Left .■ a graph with V = 8 vertices and B = 9 bonds. Right .• a ring pierced 
by a magnetic flux 6 attached to a wire (B = V = 2). 



Spectral determinant.— The spectral determinant of the Laplace operator A is formally 
defined as Sg{'~f) = det(7 — A), where 7 is the spectral parameter. This object has been 
introduced in Ref. [9] in order to study magnetization of networks of metallic wires. Despite the 
Laplace operator acts in a space of infinite dimension, Sg{'y) can be related to the determinant 
of a finite size matrix [9j : 

Note that the first product, runing over all bonds, coincides with the Dirichlet determinant 
(spectral determinant for Dirichlet conditions at all vertices). A similar decomposition was 
obtained for combinatorial graphs in Ref. [10]. The V x ^-matrix Ai has matrix elements : 

Xa + V^'^aafiCOth^/^laf.j ~ °"/^ g^^^^^ ^ " 

This expression describes the case with magnetic field : 9ap is the circulation of the vector 
potential along the wire 9ap = J^i^dxA{x). Generalization to the case of the spectral determi- 
nant of the Schrodinger (Hill) operator Sg{j) = det(7 — A + V{x)) with generalized boundary 
conditions has been obtained in Refs. [HI [12] (see also the review articles [HIE])- The result ^ 
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has been derived by two methods : (i) construction and integration of the Kernel of the operator 



(-A + 7)-i ig El El 12] : Jgdx{x\ 



\x) = 0^ In 5*5 (7). (ii) A path integral derivation [7]. 



jg uu. \ a. I I J- / — 1^7 J 

These derivations give the spectral determinant, up to a numerical factor independent on 7 (this 
is inessential for physical quantities since they are always related to d^lnSg). In the present 
article, the precise prefactor of the spectral determinant is fixed by eq. ([T]). Doing so we do not 
provide a way to determine the 7-independent prefactor from the spectrum of the graph. 

It is worth mentioning that the derivation of a (^-regularized determinant allows to define the 
prefactor of the spectral determinant within the calculation. If we denote {En} the spectrum of 
an operator O, the determinant of this latter is defined thanks to the C-function ({s) = 
as det(^ O = exp — C'(0) [E]. This approach has been used in Ref. [Tl] where result of Ref. [TTj 
for continuous boundary conditions has been obtained with a procedure fixing precisely the 
multiplicative factoiQ. 

Attachment of two graphs.— Let us consider two graphs Qi and Q2 characterized by matrices 
A^^" and generalized boundary conditions at vertices a and /3, characterized by pa- 

rameters Aq and A/3. We denote by 5*^" (7) and 82'^ (j) the corresponding spectral determinants. 







Figure 2: (a) to (c) : Attachment of two graphs Qi and Q2 (the dashed areas hide the structures 
of the graphs) : a bond {af3) is introduced, then he limit l^p ^ is taken. 

We now attach Qi and Q2 with a bond {otj3) (figure Ob). The new graph is denoted Q. The 
matrix M characterizing the new graph has the structure : 



M 








V7 







%/7 









My 



\ 



(3) 



The diagonal blocks coincide with the matrices and TWg*^ isolated graphs, provided 

a modification of the parameters describing boundary condition : A'^ = A^ -|- ^/7coth ^lap and 
\'a = \p + ^/7coth (the coth accounts for the additional wire). Then we see that 



det M = det 



det 



M'2' 



J 



sinh ^lai3 



(4) 



where J^^^ is the matrix with only one non zero matrix element equal to 1 on the diagonal 



corresponding to vertex /? : J'^i^-' = S^pS^p. Below we introduce the notation J^i = 



We connect notations of Ref. 14 with ours. Matrix A ^parameters Aa's ; det(7?(A) -\- A) ^ det ; the 



Dirichlet determinant is det(-ffi5 + ^ 



(a/3) 



2 sinh ^y^l(y.f3 



Therefore Eq. (1.1) of Ref. [IJ for the ^-regularized 



spectral determinant shows that this latter is related to eq. ((T} by Sg('y) 



Ho 



-Se(7)- 
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A' =A 

and = -^2'' that denote matrices characterizing Qi and Q2 before connection. Using 
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^ J aa \^ ^ J aa 1 + (-^1 ) aa ^/l COth ^laf3 

a httle bit of algebra gives : 

detA^ = detA4idet7W2 

x{l + ^COth^l^p[{M^')aa + {M^^)p(3] + 7(-M r' (-^2 } (6) 

In the hmit — > oo, corresponding to Dirichlet boundary condition, the determinant behaves 
Unearly with Aq,, therefore we define the spectral determinant with Dirichlet boundary condition 
at vertex a and /? as : 

Sf^'il) = lim ^tM , sY^7) = lim ^ (7) 

The Dirichlet-determinant is computed by eliminating in detTWi the column and the line cor- 
responding to vertex a. Therefore {M-i )aa = S^^^ /Si. Finally, using eqs. (I1|6I) . we obtain for 
the spectral determinant of the graph of figure [2jb : 

V7%(7)= cosh^rW [Si{l)^Sf'{^) + ^SY'\i)S2{i)] 

+ sinh^/,^ [5i(7)S2(7) + ^/75f"(7)^/752°"(7)] (8) 

At this stage it is interesting to discuss the simple case of a graph with Dirichlet boundary 
at vertices a and (3. We can take the limit A^, A^ — > oo, which corresponds to the subsi- 

tution Si AqS'P''' and 52 — > XpS^^"^ . We obtain the expected result lim;^^^ ;^^^oo = 

sinh^;,g gDir(^)gDir(^) equivalent to Spec(-A;g) = Spec(-A; gi)USpec(-A; g2)U{(^)2; n G 

W}. 

The last step of the graph attachment consists to take the limit /^^ (figure [2jc) we 
obtain : 

Sgil) = Si (7) S2°-(7) + Sf'^i^) S2(7) (9) 
which is the central result of the present article. 

Example : Ring attached to a wire (figure [ip.- If we consider for Qi a ring of perimeter L 
pierced by a flux 9 (corresponding to A{x) = we have : ^ring = 2(cosh.y7L — cos0) and 

■^ring ~ The graph Q2 is a wire of length h with Neumann boundary at its ends (for 

a vertex a of coordination number rria = 1 the case Aq, = coincides with Neumann boundary 

condition) : 5'^°*e^'^'^" ~ V^^™'^ '^wire''^''^ ~ cosh ^776 . Therefore, we recover the simple 

result : 5(7) = sinh y/^h sinh + 2 cosh y/^h (cosh ^/^L — cos 9) obtained directly from ([T]) in 
Ref. HSl. 



-1^ 



Eq. (0] is contained in eq. The result ([9]) has appeared as a limit of eq. ([8]), therefore it 

seems at first sight a particular case of this latter equation. We show now that eq. ([8]) can in 
fact be recovered from ([9]). For that purpose we proceed in two steps. First we attach a wire 
of length h to the graph Qi. The graph formed is denoted Winter and the corresponding spectral 
determinants 5inter and ^Pter' depending on the nature of the boundary condition at the end of 
the wire. Spectral determinants of the wire for the three different boundary conditions, S^^^ "^eu^ 
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•^wh-e^^"^ and 5'^°re^^''^ = "^ring' "^^re given above. Therefore, from eq. Q we obtain : 

Winter = ^1 S^Z^""'' + Sf'' S^^ = Si cosh + Sf" ^7 smh V^b (10) 

cDir c c-both Dir , c-Dir c-Neu/Dir ^ sinh p^^. 

-Winter = -^l-^wirc +"^1 -^wire = "^1 ^ ^^l COSll (H) 

In a second step we attach the graph G2 to the end of the wire of Gmter- We use again eq. Q 
from which it follows that Sg = Sinter S'^''^ + •S'inter'^'z; that precisely coincides with eq. ([8]). 

Attachment of two graphs for Schrodinger operator. Let us consider the spectral deter- 
minant for the Schrodinger operator (Hill operator) Sg{j) = det(7 — A + V{x)) with the same 
(continuous) boundary conditions as above. Let us first discuss how (jll2p are modified. Va/sixa/s), 
with Xaf3 G [0,^0/3], is the component of the scalar potential V{x) on the bond. An important 
ingredient is the solution fapixap) of the differential equation [7 — + Vai3{xai3)\fai3{xai3) = 

on the interval [0, lap], satisfying fap{0) = 1 and fapi^ap) = 0. A second independent solution of 
the differential equation is f/3a{x/3a) = fpailap — Xap) (one should not make a confusion : despite 
we use the same notation, the 2B functions fap are not the B components of a scalar function). If 

Vq/3(xq/3) = we have obviously fapix) = ^^^^ir^'^i' ^ ■ shown in Ref. [llj^ that eqs. (|l|2p 

are generalized by performing the substitution ^ coth ^1^/3 — > —f'api^) and ^.^^^^ ^ — > 
-fajjO-oLp)- The matrix M becomes Map = SapiK - Y^f^a■a^lfL^liO)] + «a/3/a/3(^a/3) e"'^"" and 
the spectral determinant takes the forrro 5*5(7) = [0(0/3) fa/si^ap)]'^ det A^. 

We consider two graphs Gi and Q2 on which lives a scalar potential V{x). If these two 
graphs are attached by a bond (a/3) where the potential vanishes [V{x) 7^ for x € ^1 U ^2 and 
V{x) = for X £ (a/3)], the structure ([3]) still holds. Therefore all results derived above are still 
valid, and in particular eqs. (j8|9p and also (|12p . 



Attachment of n graphs.— We consider a graph Q obtained by attachment of n graphs at 
the same point (figure Od). It is now easy to generalize ([9]) in order to describe this situation. 
We start from ^ : Sg = Si S^l'...^^ + Sf'' S2+...+n and use S^l'..^^ = S^'' ■ ■ ■ S^'\ Proceeding 
by recurrence, we end with : 

n 

Q \ ^ oDir oDir q oDir riDir /'1o\ 

^ ^ Dirichlet Dirichlet 



Cayley tree.— We can use (fT^ to study the case of a Cayley tree of coordination number z. We 
denote Sn the spectral determinant of a Cayley tree of depth n with Aq = Va. The spectral 
determinant for the similar graph with Dirichlet boundary at one of its end is denoted S^^"^ . We 
proceed in two steps represented on figure [3] : first we attach z — 1 such trees together, using 
(jl2p . Then we attach a wire of length b by using Q. We find : 

Sn+i = [{z - 1) Sn cosh ^b + ^S^"" sinh ^b] {S^'Y'^ (13) 
VlS^ti = [{z-l)Sn^inh^b+^S^^ co^\i^b\{S^^y-^ (14) 

^ The notations used here are shghtly different from those of Ref. [TT]. They coincide with those of Refs. |16|[5]. 

^ The fact that Sg{'y) oc det has already been demonstrated in Ref. [17) . however the remaining factor has 
been obtained in Ref. |TT] by construction of the resolvent in the graph. Note that the Dirichlet determinant 
[n(c</3) ■f'aiiQ'ap)]'^ may play a role in order to determine the full spectrum. A trivial example is the wire (with 
V{x) = 0) for which det = 7, that does not determine the spectrum. Another example is studied in detail in 
section 12 of Ref. T. The 7-dependent factor [0(11/3) /ci^(^a/3)]~^ is also important from a physical point of view 
since ^ In 5(7) (for V{x) = 0) has been shown to be related to several physical quantities [9l[7l|8]. 
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with 5i = S*^?*]^^*^" and SP" = g^^^/^^^ _ Note that in the case z = 2 the recurrence is trivially 
solved and give the spectral determinants for a wire of length nb. 




Figure 3: Cayley tree of coordination number z (here A) and depth n (here 3 before attachement 
and 4 after). 



Conclusion.— Let us come back to the initial question of the paper. Eq. ([9]) seems at first sight 
to involve only spectral information on graphs Qi, Q2 and Q, however we mentioned above that, 
in eq. ([1]), the 7- independent prefactor is not a priori fixed by spectral information. Therefore we 
can only provide here a partial answer to the initial question : given the spectral determinants 
of two metric graphs Qi and Q2, defined by ([1]), we can determine the spectrum of the graph 
Q formed by attaching the two graphs at a vertex. An interesting development would be to 
provide a relation similar to ([9|) when the spectral determinant and its 7-independent prefactor 
are constructed from the spectrum only. In the case of (^-regularization of Ref . [14] , the relation 
between (^-regularized determinant and ([T|), mentioned in a footnote above, suggests that the 
relation for (^-regularized determinants analogous to ([9]) also involves some information on the 
coordination numbers of vertices. 

The choice of continuous boundary conditions was an important hypothesis in order to 
derive eqs. (I9|12p . Another simple choice of boundary conditions, assuming continuity of the 
derivative of the field at the vertices, is examined in the appendix. This leads to a relation 
with a similar structure, eq. (I16p . A question would be to generalize the results ()12|16p to the 
case of general boundary conditions. This would require to formulate the problem with matrices 
coupling arcs jl^ since in absence of continuity of the field or its derivative, one cannot introduce 
anymore vertex variables. 

An interesting development would be to generalize ([9]) to other attachment procedures 
(graphs attached at more than one vertex). For that purpose, a helpful starting point may 
be the scattering interpretation of equation ([9]). If a graph is connected to an infinite wire, its 
spectrum is continuous and we can consider the scattering problem. A plane wave of energy 
E = —7 = k"^ sent from the infinite lead is reflected by the graph with a phase shift e''^^"''"^^'^ ^ 
given by coig[5{E)/2] = as shown in Ref. [7] (eq. (117)), where S^'''{-E) corre- 

sponds to Dirichlet boundary condition at the vertex where infinite wire is attached. We can 
associate to the two graphs Qi and Q2 two such phase shifts 8i{E) and 52{E). The spectrum of 
the graph Q obtained by attachement of Gi and G2 is given by the Bohr-Sommerfeld quantization 
condition 6i{En) + 62{En) = 2mT, that rewrites cotg[5i/2] + cotg[(52/2] = 0. Since the spectral 

^Dir 5*^^^ 

determinant vanishes on the spectrum, S{—En) = 0, this shows that S oc — I — ^ (note 
however that this argument misses a factor function of the energy ; see the footnote [3]) . The 
scattering problem has been studied for graphs with an arbitrary number of contacts (infinite 
leads) ; in particular expressions of the scattering matrix of a graph with L infinite leads is 
available in Ref. p!8] (for V{x) = 0) and [6] (for V{x) / 0). It must also be pointed that the 
question of graph attachment has been studied in Ref. [19] and in particular how to construct 
the scattering matrix of a graph in terms of subgraphs scattering matrices. All these results 
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on scattering theory in graphs might help the construction of the spectral determinant of two 
graphs attached by L > 1 vertices. 

Acknowledgments.— I thank Y. Colin de Verdiere and Alain Comtet for interesting discussions. 

Appendix : derivative continous at the vertices.— The boundary conditions discussed in 
this article ((/^(x) continuous and Oa^(^'^^(0) = Xa^a) can be interpreted as the introduction 
of a (5-potential at the vertex. They are denoted "(5-coupling" in Ref. [ITj, where "(5'-coupling" 
are also introduced. These latter correspond to continuity of the derivative : V'a/jlO) = f'a"^ 
neighbour of a and ciap fa/si^) = fJ-a V'a (the limit fia ^ oo corresponds to Neumann bound- 
ary condition ip'^ = 0). The results of the present article are easily generalized to the case of 
(5'-couplings. 

Spectral determinant.- The spectral determinant now involves the solution of the differential 
equation [7 - + Vafs{xaf3)]gat3{xaf3) = on the interval [0, satisfying g'^JO) = 1 and 

a'api^ais) = 0. The spectral determinant is given by Sg{j) = [Yl^ap) daisiL/s)]"^ detM with 
AA„/3 = ^af3[l^a " J2u O'o.'^9au{0)] - aa/3ffa/3 (^0/3) e-'^"f . Note that [ll(^ap) 9apiLfs)]~^ corresponds 
to the Neumann determinant : disconnected wires with Neumann boundary conditions fia 
00 V a). In the absence of a potential, V{x) = : 

/ 1 ^ \ Q-i^cp 
J^aP = /^a + — > ^ Oav COth ^/^lau + Qa/g ^ ■ -, ^, ■ (15) 

and Sg{-f) = {U{a(3) \/7smh^/„^) detT^. 

Graph attachment.- We consider n Graphs charaterized by spectral determinants Sk for k = 
1, ■ ■ ■ ,n. We introduce the notation S]^^^ = lim^^^oo where a is the vertex of attachment 
of the n graphs (figure [2ld). Since spectral determinants for continuous boundary conditions 
and continuous derivative have similar structures, the results obtained in this article are easily 
generalized. In particular the result (jl2p . from which other results have been derived, becomes 
for ^'-couplings 

n 

Q \ ^ cNeu cNeu a oNeu cNeu /IrN 

Dg — > O;^ ■ ■ ■ '-'fe-1 '-'A: 'Jfc+l • • • '-'n ■ UOJ 
— ..^^^^^ ^^^^^^^ ..^^^^^ 

k=\ v v 

Neumann Neumann 
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